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QUARKONIUM AS A PROBE OF 
QUARK GLUON PLASMA FORMATION 

The study of quarkonium in media has recently under- 
gone crucial developments (see e.g. [1]). Large datasets 
from heavy-ion collisions have recently become avail- 
able at RHIC displaying new features related to the quark 
gluon plasma formation characteristics like the particular 
structure of jet quenching and the very low viscosity to 
entropy ratio. In particular the quark gluon plasma looks 
more like a liquid than a plasma and the use of perturba- 
tive expansion appears to be justified only at temperature 
bigger than the deconfinement one. 

The suppression of quarkonium production in the hot 
medium remains one of the cleanest and most relevant 
probe of deconfined matter However, the use of quarko- 
nium yields as a hot-medium diagnostic tool has turned 
out to be quite challenging for several reasons. 

Quarkonium production has already been found to be 
suppressed in proton-nucleus collisions by cold-nuclear- 
matter effects, which themselves require dedicated ex- 
perimental and theoretical attention. Recombination ef- 
fects may play an additional role and thus transport prop- 
erties may become relevant to be considered. Finally, the 
heavy quark-antiquark interaction at finite temperature T 
has to be obtained from QCD [2, 1]. 

THE HEAVY QUARK INTERACTION AT 
FINITE T 

Since the publication of the famous paper [3] of Mat- 
sui and Satz arguing that color screening in a decon- 
fined QCD medium could destroy all QQ bound states at 
sufficiently high temperatures, there has been consider- 
able interest in studying the interaction between a heavy 
quark and a heavy antiquark in hot media. The expecta- 
tion was that the linearly confining potential of zero tem- 
perature would get replaced by a Debye-screened poten- 
tial at high temperatures [3]. However up to very recent 
times no proper tool to define and calculate the quarko- 



nium potential at finite T was developed. Most of the in- 
vestigations have been performed with phenomenolog- 
ical potentials whose behaviour was inspired by lattice 
calculations of the free energy. Color screening is indeed 
studied on the lattice by calculating the spatial correla- 
tion functions of a static quark and antiquark in a color- 
singlet state which propagates in Euclidean time from 
T = Oto T = 1/r, where T is the temperature (see e.g. 
[1, 4, 5] for reviews). Lattice calculations of this quan- 
tity with dynamical quarks have been also reported. The 
logarithm of the singlet correlation function is called the 
singlet free energy. In the zero-temperature Umit the sin- 
glet free energy coincides with the zero-temperature po- 
tential. Moreover at sufficiently short QQ distances, the 
singlet free energy is temperature independent and equal 
to the zero-temperature potential, while for large distance 
it shows a flattening behaviour that is thought to be re- 
lated to the screening effect. The range of interaction 
decreases with increasing temperature. For temperatures 
above the transition temperature, Tc, the heavy-quark in- 
teraction range becomes comparable to the charmonium 
radius. Based on this general observation, one would ex- 
pect that the charmonium states, as well as the excited 
bottomonium states, do not remain bound at tempera- 
tures just above the deconfinement transition, and this 
referred to as quarkonium dissociation or quarkonium 
melting. The free energy is extracted on the lattice typi- 
cally from a calculation of the quark-antiquark Polyakov 
loop correlator There are singlet and octet channels that 
are gauge dependent and one can define also an aver- 
age gauge-independent free energy. The three different 
channels show a different dependency on the QQ sepa- 
ration and thus lead to different binding energies when 
used as phenomenological potentials in the Schrodinger 
equation [5]. There is a huge literature using the singlet 
free energy or the corresponding internal energy to calcu- 
late quarkonium energies at finite T and recostructing the 
lattice meson correlation functions to understand which 
one fits better 

It is therefore very important to find a theoretical 
framework that can give us the definition of what is the 



QQ potential at finite T and a calculation tool. This is 
realized by constructing an appropriate effective field 
theory (EFT). 

For observables only sensitive to gluons and light 
quarks, a very successfull EFT called Hard Thermal 
Loop (HTL) effective theory has been derived in the past 
[6] by integrating out the hardest momenta propotional to 
T from the dynamics. However, considering also heavy 
quarkonium in the hot QCD medium, one has to consider 
in addition to the thermodynamical scales in T also the 
scales of the nonrelativistic bound state and the situation 
becomes more complicate. 

In the last few years years, there has been a remark- 
able progress in addressing the problem of outUning 
an EFT framework for quarkonium at finite tempera- 
ture and in rigorously defining the quarkonium poten- 
tial. In [11, 12], the static potential was calculated in the 
regime T^l/r> mo, where mo is the Debye mass and 
r the quark-antiquark distance, by performing an analyt- 
ical continuation of the Euchdean Wilson loop to real 
time. The calculation was done in the weak-coupling re- 
summed perturbation theory. The imaginary part of the 
gluon self energy gives an imaginary part to the static po- 
tential and hence a thermal width to the quark-antiquark 
bound state. In the same framework, the dilepton produc- 
tion rate for charmonium and bottomonium was calcu- 
lated in [13, 14]. In [15], static particles in real-time for- 
malism were considered and the potential for distances 
1/r ~ mo was derived for a hot QED plasma. The real 
part of the static potential was found to agree with the 
singlet free energy and the damping factor with the one 
found in [11]. In [16], a study of bound states in a hot 
QED plasma was performed in a non-relativistic EFT 
framework. In particular, the hydrogen atom was stud- 
ied for temperatures ranging from T <C ma^ to T ~ m, 
where the imaginary part of the potential becomes larger 
than the real part and the hydrogen ceases to exist. The 
same study has been extended to muonic hydrogen in 
[17], providing a method to estimate the effects of a finite 
charm quark mass on the dissociation temperature of bot- 
tomonium. In the next sections we report our work in the 
construction of a an EFT description of heavy quarko- 
nium at finite T. We study the real-time evolution of a 
static quark-antiquark pair in a medium of gluons and 
Ught quarks at finite temperature. For temperatures rang- 
ing from values larger to smaller than the inverse dis- 
tance of the quark and antiquark, and at short distances, 
we derive the potential between the two static sources, 
and calculate their energy and thermal decay width. We 
will see in particular that the EFT enables us to give both 
a proper definition of the heavy quarkonium potential in- 
side a hot medium and solid calculation tools to obtain it 
(at least in the weak coupUng situation). 



AN EFT FRAIMEWORK 

An EFT framework in real time and weak coupling for 
quarkonium at finite temperature was developed in [18] 
working in real time and in the regime of small coupling 
g, so that gT <^T and v ^ a^, which is expected to be 
valid for tightly bound states: T(15), J/xjf, ... . 

Quarkonium in a medium is characterized by differ- 
ent energy and momentum scales. There are the scales of 
the non-relativistic bound state: (v is the relative heavy- 
quark velocity): m, the heavy quark mass, mv, the scale 
of the typical inverse distance between the heavy quark 
and antiquark, mv^, the scale of the typical binding en- 
ergy or potential and lower energy scale. Furthermore 
there are the thermodynamical scales: the temperature T, 
the inverse of the screening length of the chromoelectric 
interactions, i.e. the Debye mass mo (~ gT in the pertur- 
bative regime) and lower scales, which we will neglect in 
the following. 

If these scales are hierarchically ordered, then we may 
expand physical observables in the ratio of such scales. If 
we separate explicitly the contributions from the differ- 
ent scales at the Lagrangian level this amounts to substi- 
tuting QCD with a hierarchy of EFTs, which are equiva- 
lent to QCD order by order in the expansion parameters. 
At zero temperature the EFTs that follow from QCD by 
integrating out the scales m and mv are called respec- 
tively Non-relativistic QCD (NRQCD) [7] and potential 
NRQCD (pNRQCD) [8? ,10]. 

We assume that the temperature is high enough that 
T ^ gT ^ mo holds but also that it is low enough for 
T and 1 /r ~ mv > mo to be satisfied, because for 
higher temperature the bound state ceases to exist. Under 
these conditions some possibihties are in order. If T 
is the next relevant scale after m, then integrating out 
T from NRQCD leads to an EFT that we may name 
NRQCDhtl^ because it contains the hard thermal loop 
(HTL) Lagrangian |6|. Subsequently integrating out the 
scale mv from NRQCDhtl leads to a thermal version 
of pNRQCD that we may call pNRQCDnxL- If the next 
relevant scale after m is mv, then integrating out mv from 
NRQCD leads to pNRQCD. If the temperature is larger 
than mv^, then the temperature may be integrated out 
from pNRQCD leading to a new version of pNRQCDnxL 
[20]. Note that, as long as the temperature is smaller than 
the scale being integrated out, the matching leading to the 
EFT may be performed putting the temperature to zero. 

The derived potential V describes the real-time evo- 
lution of a quarkonium state in a thermal medium. At 
leading order, the evolution is governed by a Schrodin- 
ger equation. In an EFT framework, the potential follows 
naturally from integrating out all contributions coming 
from modes with energy and momentum larger than the 
binding energy. For T <V the potential is simply the 
Coulomb potential. Thermal corrections affect the en- 



ergy and induce a thermal width to the quarkonium state; 
these may be relevant to describe the in medium modi- 
fications of quarkonium at low temperatures. For T >V 
the potential gets thermal contributions, which are both 
real and imaginary. 



than the Coulomb potential. The calculation has been 
performed also the in EFT framework [21] and a rela- 
tion between the Folyakov loop correlator and the sin- 
glet and octet quark-antiquark correlator [18] has been 
established in this setup. 



RESULTS OF THE EFT DESCRIPTION 

General findings in this picture are: 

• The thermal part of the potential has a real and an 
imaginary part. The imaginary part of the poten- 
tial smears out the bound state peaks of the quarko- 
nium spectral function, leading to their dissolution 
prior to the onset of Debye screening in the real 
part of the potential (see, e.g. the discussion in 
[19]). So quarkonium dissociation appears to be a 
consequence of the appearance of a thermal decay 
width rather than being due to the color screening of 
the real part of the potential; this follows from the 
observation that the thermal decay width becomes 
as large as the binding energy at a temperature at 
which color screening may not yet have set in. 

• Two mechanisms contribute to the thermal decay 
width: the imaginary part of the gluon self energy 
induced by the Landau-damping phenomenon (ex- 
isting also in QED) [11] and the quark-antiquark 
color singlet to color octet thermal break up (a new 
effect, specific of QCD) [18]. Parametrically, the 
first mechanism dominates for temperatures such 
that the Debye mass mo is larger than the binding 
energy, while the latter dominates for temperatures 
such that mo is smaller than the binding energy. 

• The obtained singlet thermal potential, V, is neither 
the color-singlet quark-antiquark free energy nor 
the internal energy. It has an imaginary part and 
may contain divergences that eventually cancel in 
physical observables [18]. 

• Temperature effects can be other than screening, 
typically they may appear as power law corrections 
or a logarithmic dependence [18, 16]. 

• The dissociation temperature goes parametrically as 

4 

TtTmelting ^ mg'i [16, 19]. 



THE FREE ENERGY IN THE EFT 
APPROACH 

In [21, 23] the Folyakov loop and the correlator of two 
Folyakov loops at finite temperature has ben calculated 
at next-to-next-to-leading order in the weak coupling 
regime and at quark-antiquark distances shorter than the 
inverse of the temperature and for Debye mass larger 



CALCULATION OF BOTTOMONIUM 
PROPERTIES AT LHC 

The EFT provides a clear definition of the potential and 
a coherent and systematical setup to calculate the masses 
and widths of the lowest quarkonium resonances at finite 
temperature. In [22] heavy quarkonium energy levels and 
decay widths in a quark-gluon plasma, below the melting 
temperature at a temperature T and screening mass w/j 
satisfying the hierarchy 



m :$> mas ^ TZT :$> ma^ ^ m^, 



(1) 



have been calculated at order ma^. This implies that 
mg^ ^ T ^ mg^ and that kT is lower than 7r7^eiting> 
i.e. quarkonium exists in the plasma. We will further 
assume that all these scales are larger than Aqcd and 
that a weak-couphng expansion is possible for all of 
them. Finally, in order to produce an expression for 
the spectrum that is accurate up to order maf , we will 
assume that [mDl{mas)Y' g- 

This situation may be relevant for bottomonium \S 
states (T(15), r\i,) at the LHC, for which it may hold 
myKiS GeV > m^tts ps 1.5 GeV > TtT w 1 GeV > 
ma^ « 0.5 GeV > mo- 

We work in the real-time formalism, that allows us to 
develop a treatment of the quarkonium in the thermal 
bath very similar to the EFT framework developed for 
zero temperature [10]. We integrate out the mass and the 
relative momentum arriving at pNRQCD. According to 
the hierarchy (1), both these scales are larger than T, 
which, therefore, may be set to zero in the matching to 
the EFT. As a consequence, the Lagrangians of NRQCD 
and pNRQCD are the same as at zero temperature. 

Integrating out T from pNRQCD modifies pNRQCD 
into hard-thermal loop (HTL) pNRQCD, pNRQCDhtl, 
1 18, 20 1 . With respect to pNRQCD, the pNRQCDnxL 
Lagrangian gets relevant modifications in two parts. 
First, the Yang-Mills Lagrangian gets an additional HTL 
part [6]. This, for instance, modifies the longitudinal 
gluon propagator in Coulomb gauge into (A:^ = k^) 



I 



(2) 



2k ko-k±iri^ 

where "+" identifies the retarded and "— " the advanced 
propagator. Second, the potentials get in addition to the 
Coulomb potential, which is the potential inherited from 
pNRQCD, a thermal part, that we call 5V. 



Potential, energy and decay width 

In the following, we wiU provide the thermal correc- 
tions to the color-singlet quark-antiquark potential, the 
thermal corrections to the spectrum and the thermal de- 
cay width, aiming at a precision of the order of ma^. 




FIGURE 1. The single line stands for a quark-antiquark 
color-singlet propagator, the double Une for a quark-antiquark 
color-octet propagator and the circle with a cross for a chromo- 
electric dipole interaction. 



Integrating out the scale T 

The relevant diagram contributing to the potential is 
shown in Fig. 1 . It reads 



3 D-1 



4-D 



(3) 



where D^y stands for the gluon propagator, hg =p^/m + 
(Xs/{6r) is the octet Hamiltonian and the loop integral 
has been regularized in dimensional regularization (D = 
4 + s and H is the subtraction point). In the loop integral, 
we integrate over the momentum region ko T and 
k'^T. Since T ^ {E — ho), we may expand 

/ / , E — h„ 



E-ho-ko + iri 
{E-hof 



+i 



i-ko+iny 



-ko + irj 
-ko + ir]Y 



{-ko + it])^ 
+ .... (4) 



The real part of the thermal correction to the color- 
singlet potential reads 



Re 5y,(r) 
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(5) 
(6) 
(7) 



where the first Une of Eq. (5) is of order g^r^T^ xE/T 
and comes from the second term in Eq.(2); the second 
one is of order g^r^T^ x {E/T)^ and come from the 



fourth term in Eq.(4), and the third one, is of order 
g2f.2j?, ijiio/T)'^, comes from a self energy insertion 
inside the gluon propagator in diagram Fig. 1; is the 
number of hght quarks. 

The imaginary part of the color-singlet potential, 
comes from the imaginary part of the self energy graph 
mentioned above and reads 



\m5V,{r)=-a,r^Tml 



2 2 
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3 
(8) 



This contribution, which may be traced back to the 
Landau-damping phenomenon, is of order g^r^T^ x 

{mo/Tf-. 

Evaluating Re 8Vs{r) and Im 5Vs{r) on a quarkonium 
state with quantum numbers n and /, we obtain the ther- 

iT) 

mal correction to the energy, SE^j, and the thermal 

(T) 

width, i', coming from the scale T: 
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where E„ 
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Integrating out the scale E 

The diagram shown in Fig. 1 also carries contributions 
coming from the energy scale E. They may be best evalu- 
ated in pNRQCDhtl by integrating over the momentum 
region ko ^ E and k ^ E and using HTL gluon prop- 
agators. Since k ^ E <^T, we may expand the Bose- 
Einstein distribution 



nB{k) = 



1 



i2r 



+ ... ; 



(11) 



moreover, since k^ E ^ mo, the HTL propagators can 
be expanded in m^/E^ <C 1. 



The momentum region k^^ E and A: ~ £ is character- 
ized by two possible momentum sub-regions. This can 
be understood by considering the integral 



J (27r)o-i Jo 2nkl-k'^-m 



"dko 



■ + ■ 



■ho-ko + iri E -ho + ko + iri 



■ (12) 



For ko ^ E and ^ ~ £, it exhibits an off-shell sub-region, 
kQ — k'^E, and a coUinear sub-region, — ^ ~ mj)/E. 
Note that, according to (1), the coUinear scale satisfies 
mg'^ ^ mj) /E ^ ing^, i.e. it is smaller than mo by a factor 
of nio/E <C 1 but still larger than the non-perturbative 
scale g^r by a factor T/E ^l. 

The thermal correction to the energy, 8E^^ , coming 
from the scale E, reads 



- — OkTmlaln^[5n^ + l-3l{l+\)\. (13) 



We note the complete cancellation of the vacuum contri- 
bution (which includes the Bethe logarithm) against the 
thermal contribution originating from the "—1/2" term 
in the expansion of the Bose-Einstein distribution (see 
Eq.(ll)). 

(E) 

The thermal width, F^/, coimng from the scale E, 
reads 



64 



32 



F^7 = 4a,^r - —aJE„ + -<r^^ 
9m 3 mn-^ao 
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where /i,o = -0.49673, hfl = 0.64070, ... .The lead- 
ing contribution is given by the first three terms, which 
are of order (X^T. This contribution to the thermal 
width is generated by the possible break up of a quark- 
antiquark color-singlel state into an unbound quark- 
antiquark color-octet state: a process that is kinematically 
allowed only in a medium. The singlet to octet break up 
is a different phenomenon with respect to the Landau 
damping. In the situation E ^ mo, the first dominates 
over the second by a factor {ma^/mD)^. 



Integrating out the scale mo 

The diagram shown in Fig. 1 also carries contributions 
coming from the energy scale mo- These contributions 
are suppressed with respect to the other terms calculated. 



TABLE 1. The pattern of IR and UV divergences 
in the quarkonium spectrum at different energy 
scales. The final result is finite. 
Scale I Vacuum I Thermal 



mOs 



mat 



I mat 



(eiR 



1 

EUV 

1 



euv 



■ —ma: 



eiR 



' ma; 



£uv 



Cancellation of divergences 

The thermal corrections to the spectrum and the ther- 
mal decay width develop divergences at the different en- 
ergy scales. These are artifacts of the scale separations 
and cancel in the final (physical) results. 

Concerning the thermal decay width, the divergence 
at the scale mal in Eq. (14), which is of ultraviolet (UV) 
origin, cancels against the infrared (IR) divergence at the 
scale T in Eq. (10). 

In the spectrum, the pattern of divergences is more 
complicated and is summarized in Tab. 1 . The table may 
be read vertically or horizontally. If read vertically, it 
shows a typical EFT cancellation mechanism: at the scale 
mOs we have non-thermal IR divergences in the poten- 
tials, these cancel against non-thermal UV divergences 
at the scale ma^ [24], the non-thermal contribution at 
the scale T is scaleless and vanishes in dimensional reg- 
ularization; thermal IR divergences at the scale T cancel 
against thermal UV divergences at the scale ma^. If read 
horizontally, it shows a cancellation mechanism that is 
famihar in thermal field theory: at the scale T, thermal IR 
divergences cancel against non-thermal IR divergences 
while non-thermal UV divergences cancel against IR di- 
vergences that appear in the potentials at the scale ma^; 
at the scale ma^ UV thermal divergences cancel against 
UV non-thermal divergences. The cancellation between 
non-thermal and thermal parts is possible because the lat- 
ter may carry temperature independent terms (see, for in- 
stance, the "-1/2" term in Eq. (11)). Note that both at 
the scales T and ma^, the spectrum is finite. 



Results 

For a quarkonium state that satisfies the hierarchy 
specified in Eq. (1) and in the following discussion, 
the complete thermal contribution to the spectrum up to 
ff{mal) is obtained by summing Eqs. (9) and (13) and 



subtracting from the latter the zero-temperature part, this 
gives 



<- c,(thermal) 27r 9 ^,2 



+ 
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3n' 



log 



-2ye 
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where L„/ are the QCD Bethe logarithms: ^1,0 — 
-81.5379, L2,o = -37.6710, . . . [25]. 

For a quarkonium state that satisfies the hierarchy 
specified in Eq. (1) and in the following discussion, the 
complete thermal width up to 0'{mal) is obtained by 
summing Eqs. (10) and (14), this gives 
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(16) 



As a qualitative summary, we observe that, at lead- 
ing order, the quarkonium masses increase quadratically 
with T, which implies the same functional increase in 
the energy of the leptons and photons produced in the 
electromagnetic decays. Electromagnetic decays occur 
at short distances ~ 1/m <C 1/T, hence the standard 
NRQCD factorization formulas hold. At leading order, 
all the temperature dependence is encoded in the wave 
function at the origin. The leading temperature correc- 
tion to it can be read from the potential and is of or- 
der n^T^/{m'^as). Hence, a quadratic dependence on the 
temperature should be observed in the frequency of pro- 
duced leptons or photons. Finally, at leading order, a 
decay width Unear with temperature is developed. The 
mechanism underlying this decay width is the color- 
singlet to color-octet thermal break-up, which implies a 
tendency of the quarkonium to decay into a continuum 
of color-octet states. 
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